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Progress in Theoretical Chemistry and Physics 

A series reporting advances in theoretical molecular and material 
sciences, including theoretical, mathematical and computational 
chemistry, physical chemistry and chemical physics 



Aim and Scope 

Science progresses by a symbiotic interaction between theory and experiment: theory is 
used to interpret experimental results and may suggest new experiments; experiment 
helps to test theoretical predictions and may lead to improved theories. Theoretical 
Chemistry (including Physical Chemistry and Chemical Physics) provides the concep- 
tual and technical background and apparatus for the rationalisation of phenomena in the 
chemical sciences. It is, therefore, a wide ranging subject, reflecting the diversity of 
molecular and related species and processes arising in chemical systems. The book 
series Progress in Theoretical Chemistry and Physics aims to report advances in 
methods and applications in this extended domain. It will comprise monographs as well 
as collections of papers on particular themes, which may arise from proceedings of 
symposia or invited papers on specific topics as well as initiatives from authors or 
translations. 

The basic theories of physics - classical mechanics and electromagnetism, relativity 
theory, quantum mechanics, statistical mechanics, quantum electrodynamics - support 
the theoretical apparatus which is used in molecular sciences. Quantum mechanics 
plays a particular role in theoretical chemistry, providing the basis for the valence 
theories which allow to interpret the structure of molecules and for the spectroscopic 
models employed in the determination of structural information from spectral patterns. 
Indeed, Quantum Chemistry often appears synonymous with Theoretical Chemistry: it 
will, therefore, constitute a major part of this book series. However, the scope of the 
series will also include other areas of theoretical chemistry, such as mathematical 
chemistry (which involves the use of algebra and topology in the analysis of molecular 
structures and reactions); molecular mechanics, molecular dynamics and chemical 
thermodynamics, which play an important role in rationalizing the geometric and 
electronic structures of molecular assemblies and polymers, clusters and crystals; 
surface, interface, solvent and solid-state effects; excited-state dynamics, reactive 
collisions, and chemical reactions. 

Recent decades have seen the emergence of a novel approach to scientific research, 
based on the exploitation of fast electronic digital computers. Computation provides a 
method of investigation which transcends the traditional division between theory and 
experiment. Computer-assisted simulation and design may afford a solution to complex 
problems which would otherwise be intractable to theoretical analysis, and may also 
provide a viable alternative to difficult or costly laboratory experiments. Though 
stemming from Theoretical Chemistry, Computational Chemistry is a field of research 
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in its own right, which can help to test theoretical predictions and may also suggest 
improved theories. 

The field of theoretical molecular sciences ranges from fundamental physical 
questions relevant to the molecular concept, through the statics and dynamics of 
isolated molecules, aggregates and materials, molecular properties and interactions, and 
the role of molecules in the biological sciences. Therefore, it involves the physical basis 
for geometric and electronic structure, states of aggregation, physical and chemical 
transformations, thermodynamic and kinetic properties, as well as unusual properties 
such as extreme flexibility or strong relativistic or quantum-field effects, extreme 
conditions such as intense radiation fields or interaction with the continuum, and the 
specificity ofbiochemical reactions. 

Theoretical chemistry has an applied branch - a part of molecular engineering, 
which involves the investigation of structure-property relationships aiming at the 
design, synthesis and application of molecules and materials endowed with specific 
functions, now in demand in such areas as molecular electronics, drug design or genetic 
engineering. Relevant properties include conductivity (normal, semi- and supra-), 
magnetism (ferro- or ferri-), optoelectronic effects (involving nonlinear response), 
photochromism and photoreactivity, radiation and thermal resistance, molecular recog- 
nition and information processing, and biological and pharmaceutical activities, as well 
as properties favouring self-assembling mechanisms and combination properties needed 
in multifunctional systems. 

Progress in Theoretical Chemistry and Physics is made at different rates in these 
various research fields. The aim of this book series is to provide timely and in-depth 
coverage of selected topics and broad-ranging yet detailed analysis of contemporary 
theories and their applications. The series will be of primary interest to those whose 
research is directly concerned with the development and application of theoretical 
approaches in the chemical sciences. It will provide up-to-date reports on theoretical 
methods for the chemist, thermodynamician or spectroscopist, the atomic, molecular or 
cluster physicist, and the biochemist or molecular biologist who wish to employ 
techniques developed in theoretical, mathematical or computational chemistry in their 
research programmes. It is also intended to provide the graduate student with a readily 
accessible documentation on various branches of theoretical chemistry, physical chem- 
istry and chemical physics. 



vt 




TABLE OF CONTENTS 



Preface ix 

Theoretical Study of Reactions Catalyzed by Acidic Zeolite 1 

X. Rozanska, R.A. van Santen and F. Hutschka 

Quantum Chemical Modeling the Location of Extraframework 29 

Metal Cations in Zeolites 

G.N. Vayssilov 

Chemical Reactions of Alkanes Catalyzed by Zeolites 39 

E. A. Furtado andM.A. Chaer Nascimento 

Ab Initio Simulations of Zeolites Reactivity 77 

J.G. Angydn, D. Parsons and Y. Jeanvoine 

Modelling of Oxide-Supported Metals 109 

M.AIfredsson, S.T. Bromley and C. R.A. Catlow 

Elementary Steps of Catalytic Processes on Metallic and 149 

Bimetallic Surfaces 

F. Illas, C. Sousa, J.R.B. Gomes, A. Clotet and J.M. Ricart 

Role of Point Deffects in the Catalytic Activation of Pd Atoms 183 

Supported on the MgO Surface 

G. Pacchioni, L. Giordano, A.M. Ferrari, S. Abbet and U. Heiz 

The Valency Effect on Reaction Pathways in Heterogeneous 199 



Catalysis : Insight from Density Functional Theory Calculations 

A. Michaelaides and P. Flu 

The Adsorption of Acetylene and Ethylene on Transition 
Metal Surfaces 

C.G.P.M. Bernardo and J.A.N.F. Gomes 

vii 



217 




TABLE OF CONTENTS 



viii 

Theoretical Approaches of the Reactivity at MgO (100) and 
Ti0 2 (110) Surfaces 

C. Minot 



Subject Index 




Preface 



Heterogeneous catalysis is a subject of outmost importance in 
the present days. Heterogeneous catalytic processes play an essential 
role in the manufacture of a wide range of products, from gasoline and 
plastics to fertilizers, which would otherwise be unobtainable or 
prohibilively expensive. Apart from manufacturing processes, 
heterogeneous catalysis is also finding other very important 
applications in the fields of pollution and environment control. 

The immense technological and economic importance of 
heterogeneous catalysis, and the inherent complexity of the catalytic 
phenonema, have stimulated theoretical and experimental studies by a 
broad spectrum of scientists, including chemists, physicists, chemical 
engineers and material scientists. 

Nowadays, due to the great sophistication achieved by many of 
the experimental techniques used in surface science studies, a great 
deal of information can be obtained, at the atomic scale, about the 
mechanism of heterogeneous catalytic processes under well controlled 
conditions. In spite of that, it is still difficult, and sometimes 
impossible, to obtain a precise picture of the mechanism of the 
reaction, at the molecular level, without any theoretical support. 

This book aims to illustrate and discuss the subject of 
heterogeneous catalysis and to show the current capabilities of the 
theoretical and computational methods for studing the various steps 
(diffusion, adsoprtion, chemical reaction, etc.) of a heterogeneous 
catalytic process. 



M.A. Chaer Nascimento 
July 2001 




THEORETICAL STUDY OF REACTIONS 
CATALYZED BY ACIDIC ZEOLITE. 



Xavier Rozanska, Rutger A. van Santen and Francois Hutschka # 

Schuit Institute of Catalysis, Technical University of Eindhoven, PO. Box 513, 5600MB 
Eindhoven, The Netherlands 

# Total FinaElf, European Research and Technical Center, B.P. 27, 76700 Harfleur, France 



Abstract After a general introduction on zeolites and their properties, we will show 

how the current quantum chemical methods allow nowadays the simulation 
of catalytic models that become more and more realistic. Predictability of 
activity, selectivity and stability based on known structures of catalysts can 
be considered the main aim of the theoretical approach to catalysis. Here 
for a particular class of heterogeneous catalysts, the acidic zeolites, and for 
a particular class of reactions, the aromatics isomerization reactions, we 
will describe a quantum mechanic study of a reaction pathway 
investigation. The cluster approach method as well as electronic structure 
periodic method will be used for this purpose. 



1. INTRODUCTION 

1.1 Zeolite Crystals 

The Greek name zeolite has been given to a class of mineral. Zeolite 
means "boiling stone" in allusion to the behavior of some silicate minerals on 
heating. This relates to a well-known property of a zeolite that is to adsorb 
large quantity of water and gases. This is permitted by the microporous 
structure of zeolite, which can accept an important amount of guest 
molecules. Zeolites are silicon oxides crystals that show various and well 
defined microporous structures (see Figure 1). It is this characteristic which 
makes them interesting molecular sieves for which application they have been 
employed for many years. 1 

The zeolite framework is constituted of the assembly of Si0 4 tetrahedra 
which link together by sharing an oxygen atom. The potential energy surface 
of the zeolite Si-O-Si angle is rather flat, 2 which explains the large variety of 
accessible structure that can be accommodated by zeolites. 3 

The zeolites have been first used as catalyst in the 1960s for alkane 
cracking reactions in petroleum industry. 3b They replaced favorably 
previously employed alumina based catalysts because of their better thermal 
and mechanical stability. Moreover, they showed higher selectivity. The 
selectivity finds its source in the zeolite micropore structure with different 

l 

M.A. Chaer Nascimento (ed.), Theoretical Aspects of Heterogeneous Catalysis, 1 - 28 . 
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consequences on the course or a reaction ' 
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Figure 1. Some examples of zeolite structures (the dots depict the Connolly surfaces). 

A first, rather obvious reason is the reactant selectivity. As the zeolite 
micropore channels have a well defined diameter, reactants bigger than this 
diameter cannot enter within the micropores to react, in case external catalytic 
active sites of zeolite have been deactivated . 43 Smaller reactants than the 
micropores will be the only ones to get involved in reactions. 

Once a reactant molecule has adsorbed within the zeolite mouth, it needs 
to diffuse toward the active sites in contact of which reactions will occur. This 
diffusion can be very dependent on the size and shape of the zeolite 
micropores as well as on the size of the reactants or products. This becomes 
especially true when the reactants or products have similar size to the 
micropores diameter . 43 After reaction, it is the turn of the products to diffuse 
away from the micropores. 

These stereo-selectivity properties do not only apply to the case of zeolite 
catalysts but also in the case of zeolite molecular sieves. 
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After reactants have adsorbed within the zeolite mouth and have diffused 
towards the actives sites, they must react. Selective adsorption depends on the 
local topology of the active site, its immediate neighborhood, and on size and 
shape of the reactant. One may say that a similar behavior as the lock and key 
principle of molecular recognition known from biochemistry applies. In this 
sense zeolites can be viewed as analog of enzymatic catalysts. 5 This 
adsorption selectivity of the reactants to the active sites favors preferential 
reaction pathways. 

Elementary reaction steps proceed via transition states. Transition state 
requires a specific geometry that differs for each reactions step. Reaction 
steps can be favored or prohibited as a function of the available space around 
the catalytic active site. 4 This effect is designated as the transition state 
selectivity. 

Whereas selective diffusion can be better investigated using classical 
dynamic or Monte Carlo simulations, 6 or experimental techniques, 7 quantum 
chemical calculations are required to analyze molecular reactivity. Quantum 
chemical dynamic simulations provide with information with a too limited 
time scale range (of the order of several hundreds of ps) 8 to be of use in 
diffusion studies which require time scale of the order of ns to s. 6 ' 7 However, 
they constitute good tools to study the behavior of reactants and products 
adsorbed in the proximity of the active site, prior to the reaction. Concerning 
reaction pathways analysis, static quantum chemistry calculations with 
molecular cluster models, allowing estimates of transition states geometries 
and properties, have been used for years. 9 The application to solids is more 
recent. 

But before continuing the discussion on the quantum chemical 
simulations, we will give a deeper description of the current understanding of 
physical chemistry properties of zeolites. We will be using experimental as 
well as theoretical studies to support the discussion. 

1.2 Zeolite Properties 

It has already been mentioned that zeolites are shape selective with 
respect to molecular adsorption. This property relates to their micropores 
structure. The zeolite framework shows a limited flexibility, which is 
essential. For instance, Yashonath et al. 10 have shown in their classical 
dynamic simulations study of molecular diffusion within zeolite micropore 
that the zeolite framework flexibility affects significantly diffusion when the 
molecules have a size comparable with the micropore size. To get an idea of 
the order of magnitude of this flexibility, one can consider the hybrid semi- 
empirical DFT periodic study of chabazite zeolite ofUgliengo et al. 11 . They 
introduced in the unit cell of chabazite Brpnsted acidic sites which are known 
to induce an increase of the volume of around 10 A 3 . 2 ' 12 This increase of the 
volume relates with the difference of volume between a SiC>4 tetraheron and a 
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A10 3 (0H) tetrahedron. Ugliengo et al. found that a linear relationship 
between the volume of the unit cell as a function of the number of introduced 
Brpnsted acidic sites is not anymore followed only above a Si/Al ratio of 3. 

In this last example, we mentioned that zeolite framework silicon atoms 
could be substituted. This substitution is naturally found in zeolites. Silicon 
atoms can be replaced with +III valence atoms, such as Af m or Ga ul , or +V 
valence atoms, such as P +v . Then, cations or anions are introduced to 
neutralize the charge of the framework (see Figure 2). The only known 
limitation of the silicon substitutions is the Lowenstein rule which states that 
substitutions cannot be found into two adjacent tetrahedra. The presence of 
cations or anions within zeolites allows the zeolite crystals to be used in ion 
exchanged processes. 14 




Figure 2. Proton bonded to bridging O-atom in the neighborhood of a silicon atom substitution 
with an aluminum atom in H-Mordenite. 

Zeolite crystals are insulators. However, they are ionic crystals, and the 
radii of the zeolite oxygen and silicon atoms is close to that of the ones 
observed for O’ and Si ' ' ions from others crystals. 15 Despite this, long range 
electrostatic contributions play a limited role, and the bonding within the 
zeolite framework is mainly covalent with the ionic bonding contributing only 
for around 10%. 16 Zeolite crystals are characterized with small dielectric 
constant: E\ is generally between 2 to 7 for full silicon zeolites. 17 Therefore, 
the short range electrostatic contributions dominate over the interaction 
energy of molecules with the zeolite host, where dispersion Van der Waals 
contribution plays the main role. As the zeolitic oxygen atoms have a large 
radius, adsorbing molecules experience interactions only with the zeolitic 
oxygen atoms (see Figure 3). 
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Figure 3. Hexane occlude in a mordenite zeolite micropore. The hydrocarbon molecule can 
only interact with zeolitic oxygen atoms. 

Kiselev et al. 19 defined a classical force field for which zeolite silicon 
atoms are not explicitly considered. The dynamic simulations they performed 
using this force field succeeded well in describing molecules diffusion within 
zeolite micropores. 19 20 

We will now consider more specifically the properties of zeolite catalysts. 

1.3 Zeolite Catalysts 

Zeolites can be used as inert support for small metallic clusters catalysts 
or can be used themselves as catalysts. For the later case, silicon substitutions 
by another cation have to occur. The catalytic activity of zeolites will be 
dependent upon the nature of the introduced cations (viz. Fe 3 , Zn 2> , Fe 2+ , 
Cu 2+ , H + , Li + , Na\ etc...) or anions (viz. F', OH' etc...). 21 Lewis basic and/or 
acid sites and/or acid Brpnsted sites are created. We will only describe the 
zeolites with silicon substituted by Al and a proton as cation to neutralize for 
the framework charge. In this case, a zeolite is a solid acid catalyst, which acts 
as more environmentally friendly acid catalyst than the classical ones. Such 
catalysts are obtained when the M + cation of a M x [Al x 'Sii. x 02] n^O zeolite is 
ion exchanged with NH4 . The heating of the NH/ exchanged zeolite induces 
NH 3 to desorb, and a proton is left behind. The proton binds to an oxygen 
atom that bridges an aluminum atom and a silicon atom (see Figure 2). As for 
other zeolitic bonds, the proton - zeolite oxygen atom shows a strong covalent 
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22 

behavior. 

23 

Solid acid zeolites, beside inducing reactions known from superacids," 
are indeed moderate acids. 24 The moderate acidity allows zeolite catalyzed 
reactions to not suffer from side chains reactions. 4 Enthalpies of reaction 
compare actually more closely with gas phase reactions than with reactions 
catalyzed by homogeneous high dipolar acids. 25 However, zeolite micropores 
have been shown to have a stabilizing effect on the carbocationic transition 
states with respect to gas phase results, which can reach 10 to 30 % of the 
activation energies. 26 This stabilization has been demonstrated to be mainly of 
short range electrostatic nature. 26,1 

A famous example, which illustrates the reaction mechanisms that happen 
with the acidic zeolite catalysts, is the chemisorption reaction of olefins (see 
Figure 4). 





secondary alkoxy: i-propoxy 

Chemisorption 

H 




primary alkoxy: n-propoxy 



Figure 4. Reaction steps involved in the chemisorption of propylene catalyzed by an acidic 

zeolite. 

This reaction has been the subject of many experimental," ’ and 
theoretical 28 studies. The reaction, which initiates from a propylene 
physisorbed to the acidic proton, leads to the formation of a more stable 
chemisorbed propylene, or alkoxy species. Such reaction has been shown 
experimentally to occur readily at room temperature within an acidic zeolite. 27 
Whereas this reaction in principle can produce two different alkoxy species 
(viz. a primary and a secondary alkoxy species), experiment reports that only 
the secondary alkoxy species can be formed. 27 This is explained by the fact 
that the formation of a transient primary carbenium ion is energetically more 
demanding than the formation of a secondary carbenium ion. 28c As already 
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mentioned, the zeolite framework stabilizes partially the carbocationic 
transition state. This stabilization is however not sufficient to allow 
carbocations to exist as stable intermediates, and carbocations become 
covalently bonded to a zeolite framework oxygen atom. 

Solid acid zeolites are used as catalysts for a large range of reactions on 
hydrocarbons. 29 They can for instance induce alkylation, transalkylation, 
isomerization and cracking reactions. Moreover, they are also used to achieve 
fine organic reactions. 5 ' 30 We will not describe all these mechanisms, and we 
will only consider the case of aromatics isomerization reactions catalyzed by 
acidic zeolite. These reactions have been used for years as a support in the 
discussion on the zeolite catalysts selectivity. 4 ' 29 They allow us to show how 
modem quantum chemical tools can be used to describe zeolite framework 
effects on the course of a reaction. These effects include the zeolite 
framework steric constraints and electrostatic contributions. Before discussing 
the reaction, we will first give a short description of the available quantum 
chemistry tools. 

1.4 Quantum Chemistry Applied to Zeolite Catalysis 

In this part, we will not give an overview of the quantum chemistry theory 
that has led to the current state of affairs, as excellent introductions can be 
found elsewhere. 815 ' 31 We will rather shortly describe which methodology can 
be followed to describe acid zeolite catalyzed reactions. 

Because of the size of the reaction centers to be considered, a 
breakthrough in quantum-chemistry has been necessary to make 
computational studies feasible on systems of catalytic interest. The Density 
Functional Theory (DFT) has provided this breakthrough. 32 Whereas in the 
Hartree-Fock based methods, mainly used before the introduction of Density 
Functional Theory, electron-exchange had to be accounted for by computation 
of integrals that contain products of four occupied orbitals, in Density 
Functional Theory these integrals are replaced by functionals that only depend 
on the electron density. This allows a consequent reduction of the 
computational costs. An exchange-correlation functional can be defined that 
accounts for exchange as well as correlation-energy. Correlation energy is the 
error made in Hartree-Fock type theories by the use of the mean-field 
approximation for electronic motion. 

The unresolved problem yet is the determination of a rigorously exact 
functional. For this reason all of the currently used functionals are 
approximate. Furthermore, one of the main problem of DFT methods is that 
Van der Waals dispersion contribution is not explicitly considered in the 
equations, which leads to some severe errors for the description of the zeolitic 
guest-host interaction. 33 We will give further details on this point later in the 
discussion. 
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1.4.1 The Cluster Approach 



We mentioned that bonding in zeolites is dominated by covalent bonding. 
Therefore, zeolite properties can be described as being mainly locally 
dependent. For instance, it is known that Brpnsted acidic sites induce 
important distortions in the zeolite framework (see Figure 3). However, these 
distortions remain local. 2 ' 34 A direct application of this property is the validity 
of the use of small fragments to describe catalytic active sites (see Figure 5). 





Figure 5. The principle of the cluster approach method. The small cluster model (right) aims to 
model an acidic zeolite catalyst (left). 

This model approximation is known as the cluster approach. 35 The 
fragment extracted from a zeolite framework is terminated with hydrogen 
atoms. Such method has been used with success for the estimate of properties 
such as frequencies of vibration, 36 NMR shielding constants, 37 or reactivity. 2 
Another mandatory reason of the use of this method was the too heavy 
computational cost to simulate larger systems. 

In recent years, progress in computer power as well as in theoretical 
methods allowed studies with increasingly larger clusters. 38 A large cluster 
model can partly allow for the zeolite framework description. Especially, 
Zygmunt et al. 2 a provided with large cluster studies which led to a deeper 
insight into reactions catalyzed by acidic zeolites. These large cluster models 
can however partly described the effect of the zeolite framework over the 
properties of interest. 39 If long range electrostatic contributions have a limited 
impact within zeolite, they still do exist, and should not be simply not 
considered as they can alter some properties. On the other hand, even small 
cluster calculations give valuable and useful information, as it will be 
described in a future section. Alternatively to the cluster approach, others 
methods have been developed. 




Theoretical Study of Reactions Catalyzed by Acidic Zeolite 



9 



1.4.2 Beyond the Cluster Approach 

A way to describe the zeolite framework at a low computational cost is to 
use quantum mechanic - molecular mechanic methods (QM/MM) (see Figure 
6). 39 With QM/MM, only the site of interest (viz. reactants and catalytic active 
site) are treated at a quantum mechanic level, whereas the zeolite framework 
is described using force field equations of molecular mechanic. 40 



QM 




MM 



Figure 6. Principle of the QM/MM method applied to a zeolitic case. Here, the adsorption of 
methanol on a Brensted site proton within a 8-membered ring zeolite. 

The quality of these methods depends on the force fields parameters, the 
way the QM and MM parts are linked, and how QM and MM parts affect 
each other. The main advantages of QM/MM are to present a limited increase 
of required computer power as a function of the size of the system. The MM 
part can be constituted by up to thousands of atoms. 12 ' 280 ' 38 A drawback is that 
it is not easy to define a priori what should be the size of the QM and MM 
parts. Ramanchandran et al. 41 observed in their periodic study that during 
transition state electron delocalization from the Br0nsted site to others zeolite 
framework oxygen atoms was an important phenomenon. Then, large QM 
part is required which makes more costly calculations. Furthermore, another 
drawback of QM/MM is the complexity of the tuning which can lead to 
misleading results. 280 

The periodic electronic structure calculations methods constitute the other 
available approach. 42 These methods require a relative high computational 
effort, beside progress in computer power allows calculations on systems of 
size of interest. 43 The advantage of periodic approach method is that the entire 
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system is described at a quantum chemical level. The disadvantages are the 
heavy computational cost which limits currently the size of the unit cell to 
systems below 300 atoms, and the usual artifacts that are associated with 
periodic boundary conditions. 44 Next, others problems originate from the 
DFT, 33 and they will be explained in a next section. 

We will give further details about the periodic electronic structure 
calculations method as we will use results obtained with this method to 
support the discussion. 45 First, we will describe how cluster approach can be 
used to investigate reactions catalyzed by acidic zeolites. Next, periodic 
electronic structure calculations will be used to enlighten the effects of the 
approximations of the cluster approach. These approximations relate mainly 
with the missing description of the zeolite framework contributions on the 
molecules involved in the reactions. Finally, by increasing the size of the 
aromatics involved in the isomerization reactions, we will show how steric 
constraints affect the course of a reaction. 

Before this, we will give details of the methods employed for the 
calculations realized in this study. 

1.4.3 Methods and Models Employed in this Study 

The cluster approach calculations have been performed using 
Gaussian98 46 with the B3LYP method. 47 This DFT method appears to be the 
optimum choice for treatment of zeolite cluster systems and presents results 
comparable with MP2 method. 43 In order to describe the Br0nsted acidic site, 
a 4 tetrahedra cluster (A[(OHSiH 3 )(OSiH 3 ) 2 (OH)) has been chosen. We 
selected the basis set d95. Geometry optimization calculations have been 
carried out to obtain local minima for reactants, adsorption complexes and 
products and to determine the saddle point for transition states (TS). 
Frequency calculations have been computed in order to check that the 
stationary points exhibit the proper number of imaginary frequencies: none 
for a minimum and one for a transition state. Zero point energy (ZPE) 
corrections have been calculated for all optimized structures. 

The Vienna Ab Initio Simulation Package (VASP) has been used to 
perform the periodic structure calculations. 423 ' 13 ' 43 All atoms have been allowed 
to relax completely within the periodic unit cell. The large 12-membered ring 
Mordenite has been used for this study as this zeolite has a relatively small 
unit cell (i.e. 146 atoms). It has previously been studied by Demuth et al. 50 . 
For our zeolite model, the Si/Al ratio is 23, and the geometry of the unit cell is 
described by a = 13.648 A, b = 13.672 A, c = 15.105 A, a = 96.792 °, J3 — 

90.003 0 and y = 90.022 °. With VASP, the energy is obtained solving the 
Kohn-Sham equation with the Perdew-Zunger exchange-correlation 
functional. 51 The results are corrected for non-locality within the generalized 
gradient approximation (GGA) with the Perdew-Wang 91 functional. 3 " This 
functional is currently the only available in VASP. This program uses plane- 
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waves and pseudopotentials, which allow a consequent reduction of 
computational costs. A cut-off of 300 eV and a Brillouin zone sampling 
restricted to the T-point have been used. A quasi-Newton forces minimization 
algorithm has been employed: convergence was assumed to be reached when 
forces were below 0.05 eV/A. The TS search method in VASP is the nudged 
elastic band (NEB) method. 53 Several images of the system are defined along 
the investigated reaction pathway. These images are optimized but only 
allowed to move perpendicularly to the hyper-tangent defined by the normal 
vector between the neighboring images. We employed up to 8 images to 
analyze transition states. When forces of the images atoms were below 0.08 
eV/A, the forces of the maximum energy image were minimized separately. 



2. TOLUENE ISOMERIZATION 

In this part, we will summarize some of our results on the investigation of 
the toluene intramolecular isomerization pathways. 45 ' 54 Both cluster approach 
and periodic approach methods have been employed which allow giving an 
illustration of the consequence of the simplistic model in the cluster approach. 
H-Mordenite (H-MOR) zeolite is used for the periodic calculations. The 
toluene molecule does not have a problem to fit within the large 12-membered 
ring channels of this zeolite. 18 Furthermore, the intramolecular transition 
states do not suffer from steric constraints. 4 It is known that intramolecular 
aromatics isomerization can proceed via two different reaction pathways (see 
Figure 7). 4 ' 55 The first route proceeds through a methyl shift isomerization, 
whereas the second route involves a dealkylation or disproportionation 
reaction which results in the formation of a methoxy species and benzene as 
intermediate. 
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Figure 7. Catalytic cycles and description of the mechanisms of the intramolecular 
isomerization reactions of toluene catalyzed by an acidic zeolite. 55 
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The initial step of this investigation is to analyze the reaction pathways 
using the low computational cost cluster approach. The small cluster model 
aims to model a zeolitic Br0nsted acidic site, and has been demonstrated to 
fill successfully this task. 2 On the other hand, a small cluster cannot describe 
the zeolite framework. By comparison of reaction pathways taking and not 
taking into account the zeolite framework, we will be able to evaluate this 
effect on reactivity. 

2.1 Cluster Approach 

The investigation of the shift isomerization reaction pathway led to the 
establishment of the reaction energy diagram plotted in Figure 8. 54 The 
geometries of the intermediates and transition states are also summarized in 
this figure. Prior to reach the shift isomerization transition state, toluene needs 
to be activated by proton attack. This activation proceeds via a transition state 
very similar to the one obtained from propylene chemisorption in acidic 
zeolite (see Figure 4). It leads to the formation of a phenoxy intermediate 
which is less stable than the physorbed toluene. This is expected from the 
chemical properties differences between propylene and toluene. 56 




Figure 8. Reaction energy diagram and description of the mechanisms of the shift 
isomerization reaction of toluene catalyzed by acidic zeolite (all data in kJ/mol). 54 

The phenoxy species is released from the cluster with no activation 
energy barrier to overcome but a constant increase in energy to a Wheland 
complex from which shift isomerization transition state takes place. With 
respect to physisorbed toluene, the activation energy to achieve this transition 
state is £ ac , = + 282 kJ/mol. In the transition state, the shifting methyl group 
occupies an intermediate position between the aromatic ring carbon atom it 
was connected to, and the carbon atom it will connect to. The shift methyl 
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group carbon atom is located in a position almost orthogonal with respect to 
the plane defined by the aromatic ring. Extra weak hydrogen bonds with the 
cluster oxygen atoms stabilize the transition state with this structure. 

For the isomerization of toluene via methoxy and benzene intermediate, 
the similar proton activation, which results in the formation of a phenoxy 
intermediate, is achieved (see Figure 9). As previously, the bond between the 
zeolitic oxygen atom and the aromatic carbon atom stretches out. A "free" 
Wheland complex is eventually reached which can reorient to favor the 
position of the toluene methyl group with the demethylation transition state. 

The geometry of this transition state can be described as a methenium 
carbocation sandwiched in between a benzene molecule and a deprotonated 
Br0nsted site. The methenium ion is planar, and its carbon atom is located 
along the line defined by the zeolitic oxygen atom to which it will be bonded 
and the aromatic carbon atom to which it was bonded. The activation energy 
which is required to achieve this reaction is £ act = + 279 kJ/mol. 




Figure 9. Reaction energy diagram and description of mechanisms of the 
disproportionation via methyl alkoxy isomerization reaction of toluene catalyzed by an acidic 

zeolite (all data in kJ/mol). 54 

This transition state results in the formation of a methoxy species and 
benzene. Benzene can desorb from the methoxy species, or can change its 
orientation. The energy level of this intermediate is + 70 kJ/mol with respect 
to physisorbed toluene. The regeneration of toluene from this intermediate 
follows the reverse transition state. 

Interestingly, the two reaction routes require similar activation energy, 
and further consideration of the entropies of activation cannot decide as well 
whether a route has preeminence on the other (the entropies of activation are 
A 298 K.Sact = - 13 and - 17 J/mol/K for the shift and via disproportionation 
isomerizations respectively). 



2.2 Periodic Approach 
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We will now describe shortly our results of the investigation of the 
toluene intramolecular isomerization catalyzed by H-Mordenite using an 
electronic periodic structure method. 45 Let us consider first the shift 
isomerization transition (see Figure 10). As it can be seen in this figure, the 
geometries of the transition states from cluster approach and from periodic 
approach are very similar. The large difference between these transition states 
is that the activation energy from periodic approach is £ act = + 179 kJ/mol 
whereas it is E m = + 282 kJ/mol from the cluster approach. Such a decrease 
of the activation energy, when a more realistic zeolite model is employed, has 
already been notified. 281 * 




Figure 10. Geometries of the shift isomerization transition states of toluene catalyzed by acidic 
zeolite as obtained from the cluster approach method (left) and the periodic structure method 

(right). 45 - 54 

This important stabilization has other consequences on the reaction 
pathway of isomerization (see Figure 11). The investigation of the activation 
of toluene by proton attack reveals a completely different picture than the one 
obtained with the cluster approach. The protonation step of toluene becomes 
an inflection point in the reaction pathway, which gives as product a 
metastable Wheland complex. 

The formation of the phenoxy intermediate turns to be unlikely to occur 
as this intermediate remains at similar energy level as observed from the 
cluster approach (i.e. + 150 kJ/mol), whereas protonation step and Wheland 
complex energy levels are around +110 kJ/mol with respect to physisorbed 
toluene. 

Zeolite framework stabilization effects uniformly all transition states and 
charged transient intermediates, and does not effect the neutral intermediates. 
Similar effect of the zeolite framework has been described by Corma et al. 58 
for another reaction. 

Exactly the same trend is obtained for the second isomerization route. An 
important consequence of this is that the activation energy barrier for the 
isomerization via methoxy and benzene transition state is found to be E m = + 
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182 kJ/mol. This means that a periodic approach method predicts as well as 
the cluster periodic method that the two alternative isomerization routes are 
competitive pathways. Qualitatively, cluster and periodic methods give the 
same result. 
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Figure II. Reaction energy diagrams of the shift isomerization reaction of toluene 
catalyzed by acid zeolite. The diagrams using black lines refer to the cluster approach results, 
and the ones using gray lines to the periodic calculations results (in kJ/mol ). 45 ' 54,57 

However, the zeolite framework effect on the reaction is not limited only 
to a stabilization of charged species. We saw already that a transition state 
from the cluster approach turns to be an inflection point when the zeolite 
framework contribution is considered. An effect exists also on transition state. 
In the case of the shift isomerization transition state, it is found an alternative 
geometry. Before protonated toluene changes its orientation with respect to 
the deprotonated Bnpnsted site, the methyl shift reaction step can be achieved 
(see Figure 12). 




Figure 12. Front view of the alternative geometry of the methyl shift isomerization transition 
state as obtained from the periodic structure calculations . 45 
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For this other transition state, the shifting methyl group is stabilized 
because of the formation of several weak hydrogen bonds with the zeolitic 
oxygen atoms located to the opposite side of the channel with respect to the 
Br0nsted site. A striking feature about this methyl shift isomerization 
transition state is that it has a similar activation energy as the other 
isomerization transition states with £ acl = + 179 kJ/mol. Of course, this result 
is strongly dependent with the local zeolite framework topology, as we will 
see in the next section. 



3. XYLENE ISOMERIZATION 

Recently, we investigated the associative alkylation reaction of toluene 
with methanol catalyzed by an acidic Mordenite (see Figures 13 and 14) by 
means of periodic ab initio calculations. 43 We observed that for this reaction 
some transition selectivity occurred, and induced sufficiently large differences 
in activation energies to explain the small changes in the para/meta/ortho 
distribution experimentally observed on large pore zeolites. 59 The para isomer 
is the more valuable product as it is an important intermediate for terphthalic 
acid, an important polymer monomer. 44 The steric constraints obtained for the 
transition state structures could be estimated from local intermediates for 
which the orientations of the toluene molecule were similar as the ones 
observed for the transition states (see Figure 14). 




Figure 13. Geometries of the associative mechanism transition states of the alkylation 
reaction of toluene with methanol catalyzed by H-MOR as obtained from the DFT periodic 
structure calculations that lead to the formation of the different xylene isomers and water. 45 
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Toluene 




Figure 14. Reaction energy diagrams of the reactions of alkylation of toluene with methanol 
catalyzed by H-MOR, which give as products p-xylene, m-xylene. or o-xylene. and water (all 
values in kJ/mol). The values correspond to energies at OK (American Chemical Society, 

200 1). 43 

We will use now the same method and Mordenite zeolite model as in the 
previous part, and investigate the isomerization of xylene isomers. 45 As 
described in the previous part, this reaction can proceed via two alternative 
routes, viz. a methyl shift isomerization, and disproportionation reactions. 
Moreover, we observed than in the case of toluene isomerization, the location 
of toluene with respect to the Brpnsted acidic site for the shift isomerization 
was of no consequence for the activation energy barrier. We will check these 
mechanisms for the three xylenes. 

It is known from experiment that xylenes do not have different adsoiption 
energies when adsorbed within fully dealuminated mordenite. 18 In our case, 
we found however that xylenes adsorption to the acidic proton show slight 
energy differences which are correlated with the local topology of the 
Brpnsted acidic site as well as the geometry of the considered xylene isomer. 
The computed adsorption energies are - 37 kJ/mol, - 30 kJ/mol, and - 33 
kJ/mol for para-xylene, meta-xylene, and ortho-xylene respectively. One 
notes that these adsorption energies are crude underestimates of the 
experimental adsorption energies, which have been reported to be around - 
130 kJ/mol. 18 This is a direct consequence of the use ofDFT method, which 
are known to not have the ability to describe properly Van der Waals 
dispersion contribution. 33 

For a study of reactivity within zeolite, it is hopefully a good 
approximation to not consider Van der Waals dispersion contribution. In a 
classical dynamic simulation of benzene and toluene within Y zeolite pores, 
Klein et al. 60 decomposed the guest-host interaction energy. They showed that 
Van der Waals dispersion contribution was almost constant along diffusion 
pathways, and that the electrostatic contributions could explain alone the 
preferred adsorption site locations of aromatics within zeolite. The importance 
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of the electrostatic contribution is enhanced when acidic sites are present. 
Furthermore, in the case ofhydrocarbon molecules reactions within zeolite, it 
is well understood that transition states are of carbocationic nature. 2 

The main reason of a proper estimate of the adsorption energy is to allow 
the comparison with experimental. Since the pioneering work of Haag, 61 it is 
well understood that in zeolites a measured activation energy is actually an 
apparent activation energy which can be expressed in the case of a first order 
reaction as: 

£.c« aPP = £ac, + (l-0>£ads (I) 

where E M is the activation energy, 6 the coverage of the reactant to the active 
site, and 2s a ds the adsorption energy of the reactant adsorbed to the active site. 

To evaluate the Van der Waals interaction energy in the case of our 
systems, we used the force field parameters provided by Deka et al., 18 and 
defined to describe the aromatics adsorption within zeolites. The Van der 
Waals interaction is then computed between the aromatic guest and the zeolite 
host, 43 ' 45 and the DFT adsorption energy is corrected with this value as: 

£ads ~ £ads + ■C'VdW (/) 

In this way, one gets the adsorption energy of the aromatic molecule adsorbed 
to the acidic proton within the zeolite micropore. It is that is used in 

(1). Interestingly, the evaluation of £vdw for all our systems gives within few 
kJ/mol the same value, and we could assume as observed from ref. 60 that the 
dispersion contribution is a constant contribution for a given adsorbate size. 
For xylene isomers, this Van der Waals correction is £ V dw = - 95 kJ/mol . 
Then, the adsorption energies for pnra-xylene, meta-xylene, and ort/zo-xylene 
are - 132 kJ/mol, - 125 kJ/mol, and - 128 kJ/mol respectively. 

Let us now consider the methyl shift isomerizations that lead from para- 
xylene to meta-xylene (see Figure 15). 




Figure 15. Front views of the transition states of shift isomerization that occur with (right) 
and without (left) reorientation of the xylene molecule after protonation as obtained from the 
periodic structure calculations (American Chemical Society, 200 1). 45 
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As it can be observed in Figure 15, in the case of the methyl shift 
isomerization transition state for which the orientation of the aromatic ring 
with respect to the Brpnsted site is as the one obtained from the cluster 
calculations, the non-participating methyl group experiences short interactions 
with the zeolitic wall. For the transition state structure that is reached without 
reorientation of the aromatic ring after it became protonated, the non- 
participating methyl group is oriented in the same direction as the large 12- 
membered ring channel, and therefore does not suffer from steric constraints 
with the zeolite wall. This results in two very different activation energies, 
which prohibit with the achievement of the transition state with an aromatic 
ring orientation as for the cluster approach transition state. The difference in 
activation energies is A£ ac , = + 63 kJ/mol, and the activation energy for the 
methyl shift isomerization without reorientation of the aromatic ring after 
protonation occurred is E aa = + 171 kJ/mol . 

In the case of the ortho to meta-xylene methyl shift isomerization 
transition states, the steric constraints are less important as the non- 
participating methyl group has more available space because of the ellipsoidal 
shape of the 12-membered ring channel (see Figure 15). Then, the activation 
energies are + 168 kJ/mol, and + 184 kJ/mol for the transition state which 
follows immediately the xylene protonation, and for the transition state which 
occurs after xylene overcame a rotation energy barrier to change its 
orientation with respect to the Brpnsted site respectively. 




Figure 16. Front and side views of the transition states and intermediate for the 
isomerization reaction via disproportionation reaction pathway of xylene molecules catalyzed 
by an acidic Mordenite as obtained from the periodic calculations (American Chemical Society, 

200 1 ). 4S 

Concerning the disproportionation isomerization reactions, the activation 
energies follow the predictions ofCorma et al. 62 using the F1SAB principle of 
Pearson. In their study, Corma et al. predicted that the activation energies for 
the alkylation of toluene with methoxy follow the ordering ortho < para. 
They estimated these data for zeolitic systems in absence of steric constraints. 
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It can be seen in Figure 16 that the non-participating methyl group for each of 
the disproportionation transition states has enough room to avoid steric 
constraints with the zeolite wall. The activation energies are + 164 kJ/mol, + 
174 kJ/mol, and + 187 kJ/mol for an alkylation/dealkylation to the ortho , 
meta, and para position respectively. 




Figure 17. Reaction energy diagrams of the intramolecular isomerization of or/Ao-xylene (top), 
para-x ylene (middle) and meta-xylene (bottom) catalyzed by H-MOR as obtained from the 
periodic structure calculations (in kJ/mol) (American Chemical Society, 200l). 4S 

These data support the experimental study of Ivanova et al. 63 or the 
theoretical cluster approach study of Blaszkowski et al. 64 . It is known that the 
alkylation of toluene with methanol can proceed via two reaction routes. The 
first one involves an associative protonation of methanol, which induces the 
methyl jump to toluene (see Figures 13 and 14), whereas the second one is a 
consecutive reaction pathway, which initiates with the formation of a methoxy 
species and a water molecule. The second step in the consecutive reaction is 
an alkylation reaction between toluene and the methoxy species. However, 
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this conclusion has to be toned down as in the present study of the alkylation 
of toluene with methoxy, the ancillary effect of water is not considered. 5 It has 
been showed that water has an ancillary and stabilizing effect on transition 
state. These can decrease the activation energy barrier for carbocationic nature 
transition state by around 20-30 kJ/mol. 5 

The different isomerization reaction pathways considered in this part are 
summarized in the reaction energy diagrams in Figure 17. 

From Figure 17, one can conclude in total that the differences in 
activation energies for the isomerization of ortho, meta, and para-xylene are 
not large enough to induce transition state selectivity. However, transition 
state selectivity exists and prevents with the isomerization to proceed through 
some reaction pathways with respect to others. This absence of transition state 
selectivity when all possible isomerization reaction routes are considered is 
indeed expected from experimental and dynamic simulation diffusion studies. 
Deka et al. 18 reported that differences in the diffusion of the three xylene 
isomers within mordenite are sufficient to explain for the distribution of 
products as observed from experimental. On the other hand, the computed 
activation energies provided with this study are in very good match with 
experimental apparent activation energies after (1) and (2) are applied. 55 ' 65 

In the next paid, we will briefly consider a larger aromatic molecule, viz. 
dimethyldibenzothiophene (DMDBT). 66 This will give us the opportunity to 
describe the zeolite selectivity that is induced by the adsorption behavior of 
reactants to the catalytic active site on the products distribution. 



4. DMDBT ISOMERIZATION 

Hydrodesulfurization (HDS) is an important process in petrochemical 
refinery as it allows decreasing the sulfur content in diesel fuel. 67 Compounds 
such as alkylated dibenzothiophene (DBT) raise some problems in HDS. The 
alkyl groups, especially when located in the 4 and 6 positions of the ring (see 
Figure 18) make alkylated DBT resistant to classical HDS catalysts because 
of the methyl groups that prevent the thiophenic sulfur atom to be in contact 
with the active site of the catalyst. 




4,6-dimethyldibenzothiophene 

46DMDBT 



3,6-dimethyldibenzothiophene 

36DMDBT 



3,7-dimethyldibenzothiophene 

37DMDBT 



Figure 18. Dimethylated dibenzothiophene derivatives. 
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However, Michaud et al. 6S and Landau et al. 69 reported that the use of an 
acid zeolite catalyst and classical HDS catalyst bifunctional catalyst was a 
very efficient way to desulfurize alkylated DBT. In this case, they observed 
that acid zeolite catalyst achieves isomerization of the more hindered sulfur 
atom alkylated DBT (viz. 46DMDBT) to compounds for which the sulfur 
atom is not anymore hindered by the methyl groups (viz. 37DMDBT). 

We will analyze the intramolecular isomerization of 46DMDBT to 
37DMDBT catalyzed by H-MOR.™ Obviously, the large size of DMDBTs 
prevents with change in the orientation of the molecule within the narrow 
pore of the zeolite catalyst (see Figure 19). Therefore, only the shift 
isomerization reactions without reorientation of the DMDBT molecules after 
they became protonated have been considered. 
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Figure 1 9. Geometries of 46, 36, and 37DMDBT adsorbed to the acidic site, and of the 
transition states of the shift isomerization of 46 to 36DMDBT and of 36 to 37DMDBT as 
obtained from the periodic structure calculations. 70 
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The geometries of the shift isomerization transition states of 46 to 
36DMDBT and 36 to 37DMDBT and of the physisorbed 46, 36, and 
37DMDBT are shown in Figure 19. In the transition states, the shifting 
methyl groups are located between the carbon atoms it was bonding/will 
bond. As for the mechanisms obtained for toluene and xylene isomers, here, 
the zeolitic oxygen atoms that belong to an 8-membered ring side-pocket of 
MOR stabilize the shifting methyl. 

The activation energy of the 46DMDBT to 36DMDBT isomerization 
with respect to adsorbed 46DMDBT is E aa = + 137 kJ/mol. It is £ ac , = + 142 
kJ/mol for the 36DMDBT to 37DMDBT reaction with respect to adsorbed 
36DMDBT. Interestingly, the activation energies for these two reactions are 
similar. The same is observed for the isomerizations from 36 to 46DMDBT 
and 37 to 36DMDBT. In this case, the computed activation energies are ■^act 
+ 155 kJ/mol and + 156 kJ/mol respectively. 

However, further consideration of the full reaction energy diagram 
reveals an interesting situation (see Figure 20). The energy levels of the 
different physisorbed DMDBTs are very different from each other because of 
more or less sterically hampered DBT methyl groups within the narrow 
mordenite pore (see Figure 19). 

TS (3-4) 6DMDBT 




AdS_46DMDBT Ads_36DMDBT Ads_37DMDBT 

Figure 20. Reaction energy diagram of the intramolecular isomerization reactions of the 
DMDBT molecules catalyzed by H-MOR as obtained from the periodic structure calculations 

(all values in kJ/mol). 70 

With respect to adsorbed 46DMDBT energy level, the energy levels of 
adsorbed 36DMDBT and adsorbed 37DMDBT are - 18 kJ/mol and - 32 
kJ/mol respectively. One can estimate using the Polanyi-Bronsted relation, 71 
which states that the changes in activation energy are proportional to the 
changes in reaction energy, that the activation energy for a methyl shift 
isomerization of DMDBT is £ ac , = + 148 kJ/mol. Then, the selective 
adsorption behavior of the different DMDBT isomers induces with transition 
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state selectivity. For 36DMDBT, it is more likely that isomerization leads to 
formation of 37DMDBT than 46DMDBT as the A£ act is 13 kJ/mol. This is the 
result that is experimentally observed. 68 ' 69 



5. CONCLUSION 

We started this paper with a rather long introduction. This was very 
important as it allows understanding zeolites minerals and their use as 
catalysts. Based on experimental information, it becomes possible to define 
inputs to the calculations. Moreover, we discussed which data can and cannot 
be accessed by means of quantum chemistry calculations. 

The study of reactivity by QM calculations concerns only a small part of a 
catalytic event, as phenomena such as macro and micro diffusion of the 
reactants and products outside and inside the zeolite micropores cannot be 
investigated. On the other hand, QM methods are the only theoretical methods 
available that provide information on the reactivity. 

The shape selectivity of zeolites in diffusion processes is now rather well 
understood. There remains a gap in the description of the reactive events. 
With methods such as the small cluster approach, the reactivity is relatively 
well described. But, such a system cannot model the zeolite framework. 
Progress in computer power permits for the emergence of increasingly larger 
cluster studies. The large cluster models provide a useful analysis of the short- 
range electrostatic effect of the zeolite framework on transition state 
structures. But it remains questionable whether the large cluster method can 
describe transition state selectivity properly. In order to keep these molecular 
clusters in shape, generally fixed geometry models are used. How far can the 
limited flexibility of the zeolite framework smooth the large energy 
differences that are due to steric constraints? It is because of this that quantum 
chemical periodic methods constitute an adequate tool to investigate selective 
reactivity within zeolites. Of course, QM/MM methods are also useful, 
especially for a more basic (MM level) analysis of the zeolite selective 
reactivity. However, considering computational cost criteria, QM/MM 
methods are often to be preferred over periodic structure QM calculations. 
Another advantage of QM/MM methods is that Van der Waals dispersion 
contributions can be modeled. 

We have shown here, with the support of selected examples, that periodic 
calculations have become possible on systems of interest. This allows for a 
deeper understanding of zeolite catalyst reactivity. The full range of the 
zeolite selectivity can now be investigated by theoretical methods. 

We have shown in the first part what are the consequences of the absence 
of the zeolite framework on a reaction. This has an important energetic effect. 
The reaction mechanisms predicted by the cluster approach appear to be very 
similar to periodic method ones. 
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Then, we investigated a reaction for which transition state selectivity is 
important. It appears that several alternative reaction pathways can be 
followed among which some lead to a minimization of the steric constraints. 
These steric constraints are strongly dependent on the local topology of the 
zeolite framework as well as on the geometry of the transition state. 

We observed a very good agreement of the computed activation energies 
with experiment. 

Finally, we presented a case for which the large size of reactants and 
products induces transition state selectivity mainly due to differences in 
adsorption. The reactivity behavior is shown to follow the Polaniy-Brpnsted 
relation. 
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1. INTRODUCTION 

Due to the presence of three-valent aluminum in tetrahedral 
framework positions instead of tetra- valent silicon, the zeolite framework of 
the alumosilicate zeolites is negatively charged and its charge is 
compensated by different cations. Since these cations are outside the lattice, 
they are easy exchangeable, which allows a wide variation of zeolite 
properties. In H-forms of zeolites, the charge-compensating cations are 
protons which are attached to framework oxygen atoms, thus forming 
bridging hydroxyl groups (Al-OH-Si). In alkali forms of zeolites as well 
other zeolites, which exchanged with metal cations, it is expected to find the 
metal ions located near A1 positions of the lattice because of the local 
negative charge there. 

Cation exchanged zeolites are successfully applied as catalysts or 
selective sorbents in separation technologies. 1 " 3 For both catalytic and 
soiption processes a concerted action of polarizing cations and basic oxygen 
atoms is important. In addition, transition metal cation embedded in zeolites 
exhibit peculiar redox properties because of the lower coordination in zeolite 
cavities compared to other supports. 4 " 7 Therefore, it is important to establish 
the strength and properties of active centers and their positions in the zeolite 
structure. Various experimental methods and simulation techniques have 
been applied to study the positions of cations in the zeolite framework and 
the interaction of the cations with guest molecules. 1 " 4 ' 8 Flere, some of the 
most recent theoretical studies of cation exchanged zeolites are summarized. 

Since the acid-basic properties as well as the local arrangement at 
individual zeolite rings depend strongly on the aluminum content of the ring, 
it is necessary to establish the actual distribution of rings with different 
numbers of aluminum atoms in a given zeolite sample. 9 This is a problem for 
zeolites with high A1 content, e.g. faujasites, while zeolites with high Si 
content contain at most one A1 center per ring. Flowever, there is still no 
certified method to measure this distribution. For zeolites with more 
complicated framework structure as MOR, MFI, FER, arises the problem 
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with various possibilities of A1 substitution in the framework positions and, 
respectively, various possible locations of the extraframework cations. 



2. ALKALI AND ALKALINE-EARTH CATIONS 



2.1 Sodium cation 

The optimized sodium cation positions in a six-ring ofFAU zeolite 
structure containing two A1 atoms in para-position (denoted as Na-Al-2p) 10 
is shown in Figure 1. This position of the cation is representative for SII 
cation position in Y and X zeolites, 11 " 13 as well as the position of Na* in Na- 
EMT zeolite. As expected, for this and the other zeolite model structures, 
Na" prefers positions near to oxygen centers bonded to A1 atoms rather than 
those of Si-O-Si bridges. Also, the cation is far from oxygen centers which 
are connected to compensating cations, an additional proton in this case. 




Figure 1. Location of a sodium cation at a model zeolite six-ring. Zeolite oxygen centers 
directed toward inside and outside the cluster are also shown. 



The mean Na-0 distances of the rings with two or three aluminum 
atoms are close to the value obtained for the simpler model Al(H)(OSiH 3 ) 3 . 6 
When Na + is present in the ring, the oxygen centers move slightly inward the 
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ring and it intereacts stronger with the oxygen centers oriented toward the 
inside the zeolite ring, denoted in Fig. 1 as O m centers. 

The calculated BE of a sodium ion at six-rings with one, two or three 
aluminum atoms are 493-519 kJ/mol. The binding energies of the cation at 
six-rings are calculated larger than those at the four-rings. For instance, the 
Na + BE at the four-ring with two aluminum atoms is by 40 kJ/mol lower 
than at the six-rings with two A1 centers. This is in agreement with the 
experimental observation that in faujasites sodium cations prefer locations at 
six-rings (SB and SI’ sites) over SHI sites, corresponding to four-ring. 11 " 13 
The binding energy of Na + at the Sill position is probably underestimated to 
some extent in our model calculations because we did not account for the 
cation interaction with neighboring four-rings. 

It is interesting to compare the relative energies of the two isomeric 
clusters with two aluminum atoms, Na-Al-2p and Na-Al-2m. Both for the 
initial negatively charged cluster and the clusters with Na + , the ring with 
aluminum atoms in para position is more stable by 21-23 kJ/mol. This 
agrees with the Dempsey's rule 14 applied to six-rings of faujasites. 

The sodium cation in the cluster with one A1 atom, Na-Al-1, is 
located almost in the plane of the ring, near the aluminum atom, at 234 and 
235 pm from the oxygen atoms of the Al-O-Si bridges. After coordination of 
probe molecules as carbon monoxide 10 or methanol, 15 the sodium ion shifts 
higher, along the axis perpendicular to the ring, 32 pm for CO and 42-47 pm 
for methanol. This could indicate that the potential energy surface for Na + 
motion perpendicularly to the ring is flat since even such a weak interaction 
as that with CO can shift the sodium ion position. After CO adsorption the 
sodium cation is above the plane of T-atoms and the BE is reduced only by 1 
kJ/mol compared to the optimized position. The energy gain from carbon 
monoxide adsorption (13 kJ/mol) easily compensates this loss of binding 
energy. Similar shift of the Na’ cation is observed also after adsorption of 
probe molecules at a cation located at the other model rings. 

The geometry optimizations of the two clusters containing two 
aluminum atoms were performed for rings whose negative charge excess is 
compensated by a proton. For both structures, Na-Al-2m and Na-Al-2p, two 
possible positions of the sodium ion were found, one at each side of the ring. 
The positions are denoted by reference to the oxygen atoms directed inward 
the ring (Oj„), syn or anti. The anti cation position is located almost in the 
plane of the oxygen atoms 0 OU i which are directed toward the outside of the 
ring; it lies 59 pm above the plane of T-atoms forNa-Al-2p (Fig. 2b) and 80 
pm above the reference plane for Na-Al-2m. Despite the fact that Na + is on 
the side of the oxygen atoms 0 OU i, it is closer to the oxygen atoms 0, n 
because they are directed inward the ring. The position of the Oj n centers 
with respect to the ring is important not only because they are near the center 
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of the ring but also because their non-bonding lone -pair orbitals are directed 
inward the ring. 




Figure 2. Location of a litium (a), sodium (b), and potassium (c) cations with respect to the 
plane of the T-atoms of the model six-ring. 



The lone -pair density of the 0 ou , centers is mainly directed toward 
the outside of the six-rings studied. Therefore, the centers O out can take part 
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in coordination of cations outside the ring or they can play a role of basic 
centers (on the walls of the supercage) in catalytic reactions or adsorption. 
Indeed, the sites with the highest proton affinity in the clusters studied 
correspond to oxygen centers, directed outward the ring. 16 

The sodium position syn is on the side of the centers Oj m at 112 pm 
below the plane of T-atoms for the cluster Na-Al-2p or 129 pm for Na-Al- 
2m. This is 60-70 pm beyond the plane of the oxygen atoms at the same side 
of the ring. In the syn configuration, Na’ is by more than 40 pm farther from 
the O oul atoms at the other side of the ring, compared to the anti position. 
The distances Na-0, n to the closer oxygen centers are by 1-5 pm shorter. The 
two sodium positions have almost the same energy for both clusters. 

The positions of a sodium cation at the cluster containing three 
aluminum atoms are similar to that with two A1 centers. The first one (anti) 
is 33 pm above the ring, the second one (syn) is 119 pm below the plane of 
T-atonis. In this ring all oxygen atoms participate in Al-O-Si bridges, but 
due to the charge compensation two of them have each a proton attached. In 
the anti position, Na + is closer to the plane of T-atoms than in clusters with 
two aluminum centers, and the cation interacts not only with two of the Oj n 
atoms (at distances of 225-227 pm), but also with one of the 0 oul oxygen 
centers (at 234 pm) which, in this cluster model, is also oriented inward the 
ring. In the syn position, Na 1 is almost at the same distance to the O in atoms, 
but much farther from the others (the mean Na-0 ou , distance increases by 40 
pm). 

Similarly to the six-ring, the four-ring with one aluminum atom has 
only two oxygen atoms connected to Al. The four-ring is smaller, Na + is 
located farther from the ring, and its interaction with both oxygen atoms of 
the Al-O-Si bridges is weaker. For this reason the binding energy of Na‘ is 
the lowest of all clusters investigated. In the four-ring with two Al centers, 
the sodium cation is at 233-235 pm from two of the oxygen atoms. The 
position is shifted in the direction opposite to the charge-compensating 
cation. 

In order to gain information on the potential energy barrier for the 
transfer of a sodium cation between the two local minima observed for three 
of the six-rings, several intermediate sodium positions at the cluster Na-Al- 
2p were partially optimized. For each structure the Na' position was fixed 
relative to the ring and the positions of the oxygen centers and the charge 
compensation proton were optimized. We found a clear an energy barrier for 
crossing the ring of only about 10 kJ/mol. The barrier occurs when Na + is in 
the plane of the oxygen atoms directed inside the ring (at about -50 pm) 
because of the short Na-0 distances. During the transfer of the cation across 
the ring, the Oj„ centers shift slightly outward to keep optimal distances. A 
similar energy check for the sodium motion perpendicular to the ring with 
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one aluminum center did not reveal other stationary points except the 
minimum described earlier. 

2.2 Other cations 

The position of the metal cation with respect to the zeolite ring 
depends strongly on the ionic radius of the cation, as is evident from x-ray 
diffraction (XRD) data. 17 The other alkali and alkaline-earth cations 18 
studied show trends similar to Na cation. They interact mainly with (i) the 
oxygen centers of Si-O-Al bridges and (ii) the oxygen centers directed 
inwards the model ring (Oj„); also, the cation was found far from the 
bridging OH group. 

The positions of the smaller cations, lithium at the Al-2p ring and 
magnesium at the Al-3 ring, are near the plane of T-atoms of the six-ring, 
shifted somewhat toward the plane of the oxygen centers O m . The shortest 
Li-O distance is 190 pm (Fig. 2a) and Mg-O distance is 203 pm. The Mg- 
O distances to the farther oxygen centers 0; n are almost the same, 220-223 
pm. The position of Mg 2 near the center of the six-ring agrees with the 
XRD data for dehydrated MgNaX zeolite (for Mg at site SII), 17 while the 
presence of water molecules results in a shift of the cation. 17 Geometry 
optimization with different starting positions of Li’ or Mg 2 ' at the six-rings 
Al-2p and Al-3, respectively, lead to the same optimized structures, 
described above, i.e. there are no syn and anti location as observed for 
sodium cation. 

The position of Ca‘ at the six-ring Al-3 is very similar to that of a 
sodium cation at the same six-ring; this has been expected from the 
similarity of the ionic radii of the two cations, 97 and 99 pm for Na and 
Ca 2 ’ , respectively. 19 The calcium cation is at 30 pm above the plane of T- 
atoms of the ring with a shortest Ca-0 distance of 229 pm. 

Two local minima were found for K at the six-ring K-Al-2p, 
similar to the syn and anti positions of a sodium cation, but farther away 
from the plane of T-atoms. The cation is located 203 pm from the plane of 
the ring when K. is at the side of the O in oxygen atoms, i.e. in syn position, 
and 171 pm from the plane, when in anti position - at the opposite side of 
the ring (Fig. 2c). At variance with Na’, both positions of K' were also 
observed at the six-ring K-Al-1, because the cation is too big to move to the 
center of the ring as does the sodium cation. The shortest K-0 distances of 
the syn position at the six-ring K-Al-2p are 260 and 278 pm, while for the 
anti position these values are slightly larger, 270 and 279 pm. 

The Lewis acidity of the sodium cation in the two positions at the 
rings with two or three A1 centers is different judged both by the Mulliken 
charge and the frequency shift of an adsorbed probe molecule carbon 
monoxide. The cation is more acidic in syn position, q(Na) = 0.61 e, and 
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Av(CO) = 38 cm while for the anti position the charge is by 0.07-0.10 e 
smaller and so is the CO frequency shift (by 6-10 cm -1 ). At variance with 
these results, the Lewis acidity of the potassium cation at both positions is 
essentially the same: the charge varies by 0.02 e and the CO frequency shift 
by 2 cirT 1 . Since Rb + and Cs + cations are farther than K from the zeolite 
ring due to their larger ionic radii , 19 one should expect even smaller 
differences between the properties of these cations in the two possible 
positions. This comparison suggests that the existence of two cationic 
positions with respect to an individual six-ring with distinguishable 
properties is a peculiarity of Na + among the alkali cations (as mentioned 
above, only one position was observed for Li + ). 



3. TRANSITION METAL CATIONS 



3.1 Copper 

Copper exchanged zeolites and especially, Cu(I)-MFI zeolite, are 
very promising catalysts for DeNOx processes in environmental protection 
and by this reason there is a great experimental 20 " 22 and theoretical 23 " 29 
interest to these materials. Sauer and coworkers 23 " 26 studies in details the 
siting, coordination, and spectral properties of Cu + ions in MFI and FER 
zeolite using combined quantum mechanics / interatomic potential function 
technique. Various possible locations of the cation have been examined first 
my molecular modeling and after that by the hybrid technique. The Cu* ion 
prefers to occupy the open space in the channel intersection and it is 
coordinated to two oxygen atoms of the A104 tetrahedron when A1 is at the 
edge of the main and sinusoidal channels. The largest binding energy of the 
cation in MFI zeolite was found for Cu + located inside a distorted six- 
membered ring (Z 6 structure, Z denotes zig-zag channels) on the wall of the 
sinusoidal channel, where it can coordinate to three or four oxygen atoms of 
the zeolite framework. The cation-oxygen distances in this case are 212 - 
227 pm . 23 The structure 12 (I denotes channel intersection) where the cation 
is coordinated to only two oxygen centers of one A 1 atom is considered to be 
the most reactive position of the cation. Especially interesting is the 
behaviour of the excited cations in singlet and triplet states 25 in different 
binding sites of Cu + . In the triplet state, the coordination of the Cu + ion to the 
zeolite framework is significantly different from that in the singlet state. The 
Cu + ion moves away from the wall of zeolite channels, and it is coordinated 
to oxygen atoms of only one AIO 4 tetrahedron. This behaviour of the excited 
states explains the appearance of two bands in the photoluminescence 
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spectra of Cu(I)-MFI zeolites. The peak with smaller emission energy are 
found for A1 located at the intersection of two channels, while the higher 
energy peak corresponds to the Cu + sites inside one of the channels. Similar 
results were observed also for Cu(I)-FER zeolite, 26 however, here the most 
stable Cu * site are located inside the zeolite channels because of smaller pore 
diameter of these materials. 

Pierloot et al. 28 considered coordination of Cu 2 ion at cluster models 
representing six-ring sites with different Al contents. As expected, the 
optimized structures indicated a strong preference of the Cu 2 ' ion for 
coordination to oxygen centers connected to Al rather than Si. They have a 
planar - four-fold oxygen coordination in the six-rings. Depending on the 
number and relative positions of the aluminum atoms in the ring, two distinct 
coordination modes were distinguished with respect to the positions of the 
aluminium atoms in the ring. The optimized structures were used in 
clarification the peculiarities of the UV and ESR spectra of Cu 2+ -A, Cu 2+ -Y 
and Cu 2+ -ZK4 zeolites. 

3.2 Zinc 

The Zn-exchanged zeolites was studied theoretically by van Santen 
and coworkers using different zeolite models. 30 " 32 The Zn 2+ cations were 
located at four-, five-and six-membered rings, as well as on two four-rings 
with a common oxygen bridge. The binding energy of the cation is higher to 
the six-ring and lower to the four-ring. The Zn 2+ cation is the most exposed 
to probe molecules when situated at the four-ring of zeolites. A structure of 
the type [Zn-O-Zn ] 2 was found to be one of the stable forms of the Zn 
cation on zeolites with a low Si/Al ratio. This site is more 
thermodynamically stable than Zn 2 at the five-ring, which suggests that [Zn- 
O-Zn]* species likely exist in MFI zeolites. Owing to its size, this oxide 
cluster may only be stable in special positions inside the pores. 

33 Rhodium 

Coordination of Rh + in FAU zeolites was modeled by Goellner et 
al 33 Three types of cationic sites typical for faujasites were considered - a 
four-ring, a six-ring, or a three-hollow position next to the Al center, denoted 
as T5 and modeled as a four-membered ring with an additional OSiFh group 
attached at the Al atom. The calculations show that at four-rings and at six- 
rings, a Rh* ion is bonded to the two oxygen centers of the ring connected to 
the Al center, as a consequence of their high basicity; 16 the calculated Rh-0 
distances are 220-225 pm. At the four-ring, Rh is at a distance of 220 pm 
from the plane ofT-atoms and interacts with two of the oxygen centers that 
are on the same side of the ring. At the six-ring the cation is exactly in the 




